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This study of elastic and plastic deformation of graphene, silicene and boron nitride (BN) honey- 
comb nanoribbons under uniaxial tension determines their elastic constants and reveals interesting 
features. In the course of stretching in the elastic range, the electronic and magnetic properties can 
be strongly modified. In particular, it is shown that the band gap of a specific armchair nanoribbon 
is closed under strain and highest valance and lowest conduction bands are linearized. This way, 
the massless Dirac fermion behavior can be attained even in a semiconducting nanoribbon. Un- 
der plastic deformation, the honeycomb structure changes irreversibly and offers a number of new 
structures and functionalities. Cage like structures, even suspended atomic chains can be derived 
between two honeycomb flakes. Present work elaborates on the recent experiments by Jin et al, 
Phys. Rev. Lett. 102, 205501 (2009) deriving carbon chains from graphene. Furthermore, the 
similar formations of atomic chains from BN and Si nanoribbons are predicted. 

PACS numbers: 62.25.-g, 61.48.De, 61.46.-w, 73.63.-b 



I. INTRODUCTION 

For last two decades honeycomb structured materi- 
als have dominated nanoscience. The unique orbital 
symmetry of the honeycomb structure underlies several 
exceptional properties of carbon based nanomaterials, 
such as fullerenes, nanotubes, graphene and its quasi 
one dimensional ribbons. While 7r-orbitals are respon- 
sible for th e un usual electronic and magnetic properties 
of graphene^l its planar flexibility but high in-plane 
strength is achieved by sp 2 -hybrid orbitals. For exam- 
ple, charge carriers in graphene behave like a massless 
Dirac fermioniP due to electron and hole bands showing 
linear crossing at the Fermi level. These bands are de- 
rived from 7T- and 7r*-states. An unpaired 7r-state leads 
to a local magnetic moment in a nonmagnetic honey- 
comb structure^. Also the strong overlap between near- 
est neighbor 7r-orbitals assures the planar stability of 
graphene and BN. Silicene, a honeycomb structure of Si, 
lacking such an overlap is stabilized only by puckering^. 

The recent spectroscopy techniques to identify and 
quantify the strain profiles in graphene have shown 
that the Raman peaks shift considerably under the in- 
plane strairP providing a fundamental tool for graphene- 
based micro/nano mechanical systems. Some theoretical 
studies^! have also shown that local and uniform strain 
can be an effective ways of tuning the electronic structure 
and transport characteristics of graphene devices to gen- 
erate confined states, quantum wires, and collimation. A 
related step was given by Kim and collaborators^ who 
have developed a simple method to grow and transfer 
large scale, high-quality stretchable graphene films us- 
ing CVD on nickel layers, which might enable numerous 
applications including in large-scale flexible, stretchable, 
foldable transparent electronics. 

This study investigates the stretching of quasi one- 
dimensional nanoribbons (NRs) of graphene, silicene and 



BN and predicts that they attain additional functionali- 
ties by changing to a number of structures with interest- 
ing electronic and magnetic properties. While the elastic 
deformation with harmonic and anharmonic ranges and 
sudden yielding points are common to all NRs, the ab- 
sence of sequential elastic deformation stages ending with 
order-disorder structural transformation each leading to 
stepwise neckin g con stitutes their prime difference from 
metal nanowire d 10 * 11 !. Cage structures of large polygons 
can be generated. In particular, atomic chains of C, BN 
and Si between honeycomb flakes can form under cer- 
tain circumstances. The synthesis of these ultimate one- 
dimensional atomic chains, which can form various stable 
geometries^, are expected to be an essential step towards 
future "active" nanotechnology applications. We believe 
that our predictions are relevant for the current studies 
investigating the properties of strained graphene^. 



II. MODEL AND METHODOLOGY 

We have performed first-principles plane wave calcu- 
lations within Density Functional Theory (DFT) using 
PAW potentials!^ The exchange correlation potential 
has been approximated by Generalized Gradient Approx- 
imation (GGA) using PW9lG^l functional both for spin- 
polarized and spin-unpolarized cases. Recently, spin- 
polarized calculations within DFT have been carried 
out successfully to investigate magnetic properties of va- 
cancy defects in 2D honeycomb structures. Also inter- 
esting spintronic properties of nanoribbons have been re- 
vealed using spin-polarized DFT 15 . The success of spin- 
polarized DFT calculations has been discussed in Refill 

All structures have been treated within supercell ge- 
ometry using the periodic boundary conditions. A plane- 
wave basis set with kinetic energy cutoff of 400 cV has 
been used. The interaction between monolayers in adja- 
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FIG. 1: (Color online) (a) Two dimensional honeycomb structure with primitive unit cell, it— and a— bonds are schematically 
described, (b) An armchair graphene nanoribbon AGNR(9) has mirror symmetry with respect to its axis and Na—9 dimer 
bonds across the unit cell specifying its width. The unit cell is delineated by a rectangle with dotted edges. The supercell consists 
of five unit cells with Co = 5 x a. (c) Same for AGNR(IO), which lacks the mirror symmetry, (d) Response of AGNR(9) to 
uniaxial tension examined in the supercell; variation of the strain energy Es and tension force Ft with shaded region indicating 
the plastic range; atomic structure for e =0.20 just before breaking. Isosurfaces for the difference of spin up (blue/dark) and 
spin down (yellow/light) charge densities, Ap — p(t) — calculated for broken pieces show antiferromagnetic order of zigzag 
edge states. (c) Same for AGNR(IO). Monatomic carbon chains connect magnetic pieces, (d) Armchair boron-nitride NR. (e) 
Armchair silicene NR. 



cent supercells is examined as a function of their spac- 
ing. Since the total energy per cell has changed less than 
1 meV upon increasing the spacing from 10 A to 15 
A , we used the spacing of ~ 10 A in the calculations. 
In the self-consistent potential and total energy calcula- 
tions the Brillouin zone (BZ) is sampled in k-space us- 
ing Monkhorst-Pack scheme by 25x1x1 for nanoribbons. 
This sampling is scaled according to the size of superlat- 
tices. All atomic positions and lattice constants are op- 
timized by using the conjugate gradient method, where 
the total energy and atomic forces are minimized. The 
convergence for energy is chosen as 10 -5 eV between two 
steps. Numerical plane wave calculations are performed 
by using VASP packageP^ESl 



III. ELASTIC AND PLASTIC DEFORMATION 
OF NANORIBBONS 

Nanomechanics of both armchair and zigzag NRs of 
graphene, silicene and BN is explored by calculating the 
mechanical properties as a response to stretching along 
the axis of the ribbon. Mechanical properties are re- 
vealed from the strain energy Es = £r(e) — £r(e = 0); 
namely, the total energy at a given axial strain e minus 
the total energy at zero strain. Segments of quasi-one 
dimensional NRs are treated within supercell geometry 



using periodic boundary conditions. Each supercell con- 
tains n unit cells of the ribbon and hence has the lattice 
constant along the axis of ribbon, c = na, a being lat- 
tice constant for the primitive unit cell of NR. The strain, 
e = Ac/ Co, corresponds to a stretching, where the lattice 
constant of the strained supercell equals c = Co+Ac. The 
stretching of the ribbon is achieved by first increasing 
the optimized lattice constant in increments of Ae =0.01 
(namely c — > c + Ae x cq) and by uniformly expanding 
the atomic structure obtained from previous optimiza- 
tion. Subsequently the atomic structure is reoptimized 
keeping the increased lattice constant c fixed and the 
corresponding strain energy is calculated. This process 
is repeated after each increment of Ae. Then the ten- 
sion force, Ft = —dEs(e)/dc and the force constant 
k = d 2 Es/dc 2 are obtained from the strain energy. Ow- 
ing to ambiguities in denning the Young's modulus of 
2D honeycomb structures , one can use in-plane stiffness, 
C = (1/Aq) ■ (d 2 E s /de 2 f^2l. Here A is the equilib- 
rium area of the supercell. The in-plane stiffness can be 
deduced from k by defining an effective width for the 
ribbon. 

The variation of strain energy, tensile force and the 
corresponding atomic structure of selected armchair NRs 
are presented in Fig. [T] as a function of e. Here the bonds 
parallel to the direction of applied tension is stretched 
more than those in other directions. As a result, the 



3 




1.30 A 



FIG. 2: (Color online) Stretching of a finite size segment of the 
bare armchair graphene NR (Na=W) between two tapered 
ends, (a) The atomic structure of the NR corresponding to 
e=0 and e=0.17. Note that a tear starts at one edge of the 
stretched NR. (b) In the course of stretching, a hexagon con- 
nected to the chain is transformed to a pentagon by yielding 
a single carbon atom to the chain, which, in turn, becomes 
longer. In the left panel, one atom is incorporated in the chain 
when one bond of the hexagon is broken (c) A suspended 
chain comprising 12 carbon atoms is derived from graphene 
NR at e—0.40. The difference charge density of different spins 
states, Ap(f , \.) is also shown, (d) The chain structure with 
alternating short and long bonds are highlighted. The total 
charge density pr of both spins is denser around short C-C 
bonds. 



hexagonal symmetry is disturbed, but overall honeycomb 
like structure is maintained. The elastic deformation 
is reversible and stretched ribbons can return to their 
original geometry when the tension is released. In the 
harmonic range, the force constant is calculated to be 
re=176, 30 and 144 N/m, for armchair graphene, silicone 
and BN NRs having A^=10, respectively. Similarly, the 
calculated in-plane stiffness for the same ribbons are, 
respectively, C=292, 51 and 239 N/M. Notably, k and 
C values of silicene are lowest due to sp 2 orbitals de- 
hybridized as a result of puckering. Cohesive energies 
of graphene, silicene and BN nanoribbons with Na=10 
(namely 19.77, 10.76, 18.80 eV per atom pair, respec- 
tively) show trends similar to the corresponding in-plane 
stiffness values. We also calculated the in-plane stiff- 
ness values of 2D graphene, silicene and BN honeycomb 



structures to be, respectively, C=335 N/m (the reported 
experimental value^for graphene: C=340±50), 62, 258 
N/m. Because of edge effects of NRs, their stiffness val- 
ues are relatively smaller than those of 2D structures. 

The elastic range ends at the yielding point with the 
corresponding critical strain, , where the strain energy 
drops suddenly. It should be noted that ey and plastic 
deformation of nanoribbons are expected to depend on 
the ambient temperature, unit cell size, crystalline de- 
fects such as vacancy and time rate of change of stretch- 
ing. The stochastic nature of deformation is avoided to 
some extent by carrying out a slow or " adiabatic" stretch- 
ing as explained above, whereby the ribbon is stretched 
in small increments and the structure is optimized af- 
ter each increments in elastic and plastic ranges. The 
stretching of nanoribbon is carried out at 0° K for three 
different values of Ae, namely 0.05, 0.01, 0.002. The 
value of ey did not changed for Ae=0.01 and 0.002. We 
therefore concluded that already Ae=0.01 corresponds 
to a very slow or "adiabatic" stretching and reveals the 
bare response of nanoribbon without the effect of time 
rate of change of stretching. The effect of temperature is 
investigated by performing ab-initio molecular dynamic 
calculations (lasting 2 ps with time steps of 2xl0~ 15 sec- 
onds) for Na=W. Results indicate that ey decreases 
with increasing temperature. Hence, owing to the soft- 
ening of acoustical phonons, ey=0.22 corresponding to 
T=0 K is reduced to 0.16 at T=600 K. We also found 
that the results are converged if the supercell size n > 5. 
The presence of a vacancy defect in the ribbon speeds up 
the yielding by decreasing the value of ey . 

The plastic deformation stage following the yielding 
point is the crucial part in the stretching of NRs having 
honeycomb structure and hence is the focus of this work. 
The response of the ribbon to the strain after the yielding 
point is material and geometry specific. Having deter- 
mined various effects, which possibly change the value of 
ey, we examine the structure in the range of plastic (ir- 
reversible) deformation through "adiabatic" stretching. 
The armchair graphene NR with N A =9, i.e. AGNR(9) 
has a mirror symmetry relative to its x-axis. This NR is 
broken into graphene patches having equilibrium honey- 
comb structure just after ey =0.21 in Fig.[T|a). Whereas 
the behavior of AGNR(10) (which lacks the mirror sym- 
metry) is dramatically different. The ribbon is torn into 
two pieces (patches), which are connected by an atomic 
chain. In the plastic range, the strain energy E$ increases 
slightly with strain, but eventually decreases each time 
when a sp 2 C-C bond of the zigzag edge of the patch 
is broken and a C atom is incorporated into the chain 
from the graphene patch. This important result actu- 
ally predicts the recent finding by Iijima and his col- 
laborators, who derived monatomic carbon chain from 
graphene^]. Carbon atomic chains identified as cumu- 
lene (having double bonds) or polyyne (consisting of 
altern ating tr iple and single bonds) have been studied 
earlier^ 23 ! 24 !. The chain structure was found to be sta- 
ble and linear owing to the strong overlap of 7r-orbitals 
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between adjacent atoms. While infinite chain is subject 
to a Peierls distortion, bond alternation and bond length 
variation depends on the number of carbon atoms in a 
finite carbon chain. The character of the covalent sp+ir- 
bonding between carbon atoms underlies their unusual 
chemical, mechanical and quantum transport properties. 

Our results show that the tearing of the armchair 
graphene NR and hence the formation of a carbon chain 
is promoted by a vacancy defect. Even more surpris- 
ing is that not only graphene, but also puckered silicene 
and flat BN armchair ribbons are plastically deformed 
to form atomic chains between patches. Since BN hon- 
eycomb structure is already synthesized, the present re- 
sults concerning the formation of BN chains is impor- 
tant and hence is yet to be realized experimentally. In 
contrast to graphene, pieces (patches) torn from the BN 
and silicene ribbons allow also different types of poly- 
gons ranging from trigon to heptagon. In particular, a 
large hole in the silicene pieces is reminiscent of the cage 
structure as if a 2D analogue of metal-organic frameworks 
(MOFs). We also note that nonmagnetic armchair NRs 
attain spin-polarized (magnetic) ground state after they 
are broken into small pieces having zigzag edgeiPl This is 
demonstrated by isosurfaces for the difference of spin up 
and spin down charge densities, Ap(|, J,) = p(f) — p(|) 
inFig.[l] 

The behavior of a finite segment of bare armchair 
graphene NR under uniaxial tension between its two ends 
is also presented in Fig.[2ja). The tear, which initiates at 
one edge, propagates until the other edge and eventually 
the chain formation sets in. Usually triangles of atoms 
are formed at the region of junction of the nanoribbon 
pieces and the atomic chain. Upon stretching, the apex 
atom of the triangle is incorporated in the chain leav- 
ing behind a pentagonal (or broken hexagonal) ring as 
shown in Fig. j2^b) . This way, a carbon atomic chain 
with alternating long and short bonds is suspended be- 
tween two graphene pieces and grows by sequential im- 
plementation of atoms from these pieces to the chain. 
Here we point out an important difference between planar 
honeycomb NRs and metal nanorods both stretching in 
the plastic range. Experimental studies®' and theoretical 
simulations^ have demonstrated that metallic nanorods 
elongate in terms of sequential elastic and yielding stages. 
At the end of each elastic stage, the nanorod undergoes 
an order-disorder transformation through the yielding 
stage resulting in a necking. Thus the onset of a sub- 
sequent elastic stage progresses through a smaller cross 
section. These structural changes were revealed by in- 
situ measurements of two terminal conductance 10 . As for 
NRs here, the order-disorder transformation occurring at 
the end of each elastic stage and leading to necking is ab- 
sent. Necking of a NR takes place through tearing. The 
transmission coefficient, T of carbon chain, which is cal- 
culated self-consistently using non-equilibrium Green's 
Function method reflects the combined electronic struc- 
ture of central region and connecting electrodes. Results 
to be published elsewhere indicate that the implementa- 



tion of a single carbon atom from hexagon to the chain 
is reflected to the variation of T under constant bias. 



IV. VARIATION OF ELECTRONIC AND 
MAGNETIC PROPERTIES WITH STRAIN 

After massive structural changes taking place for e > 
ey, a NR attains new properties which are absent in its 
equilibrium state. For example, the ribbon achieves a 
higher chemical reactivity, because of the unsaturated 
bonds protruding from atoms having lower coordination. 
Cumulene by itself is very reactive. Not only mechanical 
properties and atomic configuration, but also the elec- 
tronic and magnetic properties of nanoribbons can be 
modified through stretching as illustrated in Fig. [3] De- 
pending on their widths, symmetries and materials the 
band gaps of nanoribbons exhibit significant variations in 
the elastic deformation range, but usually they vanish in 
the plastic range. For example, variations of the band 
gaps of hydrogen saturated AGNR(9) and AGNR(IO) 
with e in the elastic range are rather different. The band 
gap of AGNR(9) is ~0.7 eV at e=0.0, but it increases 
with increasing strain up to e=0.07, but passes through 
a maximum and subsequently decreases to vanish at ey 
as seen from Fig. |3ja). The band structure near the 
band gap and isosurfaces of charge densities of the low- 
est (highest) two states in the conduction (valence) band 
are presented in Fig. [3^b) for e = 0.0. According to the 
electronic energy bands calculated as a function strain, 
while the highest valence band is lowered, the lowest con- 
duction band is raised with increasing e This, normally, 
increases the band gap. On the contrary, second highest 
(lowest) valence (conduction) band is raised (lowered) 
with increasing strain e. Consequently, the band gap of 
AGNR(9) first increases up to e=0.07, but decreases for 
e >0.07, where the second bands cross the first ones and 
dip in the gap. At the end the orderings of the first and 
second valence and conduction bands are switched. In 
Fig. |3jc)-(d), the character and orbital compositions of 
the first and second valence and conduction bands are 
reversed in AGNR(IO). Consequently, the variation of 
band gap with the strain is reversed in AGNR(IO) rel- 
ative to AGNR(9). While the highest (lowest) valence 
(conduction) band is raised (lowered), second valence 
(conduction) band is lowered (raised) and their disper- 
sions is decreased with increasing e. This way the band 
gap decreases and is eventually closed within 3 meV at 
e ~ 0.09. Even more remarkable is that the first conduc- 
tion and valence bands, which are closed at k— ► 0, are 
linearized for 0. This means that the band energy has 
linear dispersion even for a small bias AU; namely, 

E(k) - (E F ± AV) = Ck, (1) 

where C is a constant. Accordingly, hydrogen termi- 
nated armchair nanoribbons AGNR(10), which is nor- 
mally a semiconductor, behave like a 2D graphene and 
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FIG. 3: (Color online) Variation of the energy band gaps of AGNR(9) and AGNR(10) with the strain from e=0.0 to e=0.25. (a) 
The band gap of AGNR(9) first increases with increasing strain in the elastic range, passes through a maximum, then decreases 
and eventually vanishes after the yielding point, (b) The band structures for different strain values and the isosurfaces of charge 
densities of lowest two conduction and highest two valance bands for zero-strain configuration, (c) Variation of the band gap 
of AGNR(10) with strain displays a reverse trend relative to AGNR(9). (d) same as (b) for AGNR(10). 



hence have carriers, i.e. holes or electrons, with mass- 
less Dirac fermion character at a specific value of e. This 
result is somehow unexpected but has important con- 
sequences: The massless Dirac fermion behavior in 2D 
graphene originating from linear band crossing at K- 
points at the Fermi level disappears when its size is finite. 
For example, massless Dirac fermion behavior is absent 
in an armchair nanoribbon, which is a nonmagnetic semi- 
conductor at e=0.0. Whereas the massless Dirac fermion 
behavior would be highly desirable to achieve high carrier 
mobility in these nanoribbons. Realization of linearized 
bands as predicted in Fig. [3^b) is encouraging for vari- 
ous electronic applications, especially for fast nanoelec- 
tronics. We also note that different response of different 
bands to the tensile strain gives rise to a metal-insulator 
transition. For e > 0.09, where the character and orbital 
composition of the highest valence and lowest conduction 
band switches the band gap starts to open and to increase 
with increasing tensile strain, but vanishes suddenly at 
e ~ e Y 

Because of zone folding the character and the orbital 
composition of the bands at the edges of valance and 
conduction band are switched by going from AGNR(9) 



to AGNR(IO). Therefore, similar effects can be observed 
in other families of AGNR having different widths, Na- 
On the other hand, different response of different bands 
to the strain is an interplay between bond length and 
bonding (antibonding) bond energies. In fact, we dis- 
tinguish two different bonds; namely lateral (parallel to 
the ribbon axis) and tilted bonds, which have different re- 
sponse to e. As e in the elastic range increases, the length 
of the tilted bonds in AGNR(9) first increase, passes a 
maximum and decreases reminiscent of the variation of 
the band gap in Fig. ^a) . The values of the strain corre- 
sponding to maximum of band gap and maximum length 
of tilted bonds are close. The elongation of the tilted 
bonds differ depending on their location relative to the 
edge. As for the the lengths of the lateral bonds, they 
increase with increasing e, even if their elongations differ 
depending on their positions relative to the edge of the 
ribbon. 

Fig. [4] shows the variation of the band gaps and mag- 
netic properties for other nanoribbons. In contrast to 
AGNR(IO), the wide band gap of ABNR(IO) in Fig.ga) 
doesn't change considerably for < e < 0.20 in the elas- 
tic range. However, the wide band gap of ZBNR(9) de- 
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FIG. 4: (Color online) (a) Variation of the band gap of the 
armchair BN nanoribbon, ABNR(IO), with strain e. (b) 
Zigzag BN nanoribbon ZBNR(IO). (c) The magnetic moment 
of AGNR(IO) having a single vacancy, [i jumps from 1 /j,b to 
8 fin after the NR is torn suddenly from one edge, (d) An- 
tiferromagnetic ZBNR(9) attains spin-polarized ground state 
after the yielding point, (e) Magnetic moment of a single 
edge atom of ZGNR(9) (which is antiferromagnetic in equi- 
librium) increases with strain in the elastic range and then 
falls suddenly at the yielding point. 



creases steadily from 4 to 0.3 eV between < e < 0.20 
as shown in Fig. 4jb). As a result of plastic deformation 
either a spin polarized state is induced or the existing 
magnetic state is modified. The magnetic moment of 1 
/Us of AGNR(10) having a single carbon vacancy in a 
supercell of 5 unit cells increases to 8(1b in the plastic 
range as shown in Fig. |4^c). This dramatic increase of 
the magnetic moment is attained by the severe modifica- 
tion of the honeycomb structure after the yielding point. 
FigurepJd)shows that the antiferromangetic spin state of 
ZBNR(9) changes into a ferromagnetic state. The mag- 
netization of edge atoms of zigzag graphene nanoribbon 
ZGNR(9) increases in the elastic range, but falls to a 
lower value as shown Fig. |4]Je). These results suggest 
that uniaxial strain can be used to monitor the electronic 



and magnetic properties of ID nanoribbons both in the 
elastic and plastic deformation ranges. 



V. CONCLUSIONS 

In this study the elastic constants of graphene, sil- 
icene and BN zigzag and armchair nanoribbons are de- 
termined their unusual featurs are revealed under tensile 
stress. Their atomic, electronic and magnetic structures 
are examined under elastic and plastic deformation range 
attained by "adiabatic" stretching. We found that in 
the course of elastic stretching the electronic structure 
of these nanoribbons are strongly modified. The varia- 
tion of band gap is sample and materials specific. For 
example, we showed that the variations of the band gaps 
of AGNR(9) and AGNR(10) with strain display reverse 
trends. The variation of band gap involves a complex in- 
terplay of zone folding, diverse elongation of lateral and 
tilted C-C bonds, and different orbital composition of 
the first and second valance and conduction bands. In 
particular, the closing of gap and linearization of highest 
and lowest conduction bands of a hydrogen saturated, 
armchair nanoribbon may have important implications 
, as such that massless Dirac fermion character can be 
realized even in semiconducting armchair nanoribbons. 
Unusual responses of band gaps to the strain are also 
obtained in different types nanoribbons. 

The ending of elastic range and the onset of plastic de- 
formation leading to diverse structural deformations and 
magnetic states in periodic and finite size nanoribbons 
is another interesting outcome of this study. Structures 
having large holes are reminiscent of a metal-organic 
frameworks, MOFs. We showed that long monatomic 
carbon chains can form in the course of stretching. Our 
prediction, that suspended atomic chains can also be de- 
rived from BN and silicene nanoribbons under stretching 
in the plastic range is yet to be realized experimentally. 
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